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1. A very short introduction into the history of rogue waves. 

The goal of this paper is to show that the Kadomtsev-Petviashvili equation provides a 
rough description and classification of large waves and predicts the existence of so-called 
" rogue waves" . 

Rogue waves (also known as monster waves, freak waves, killer waves, extreme waves, 
abnormal waves, etc.) are large ocean surface waves that occur far out in sea and appear 
spontaneously supposedly out of nowhere, even in calm seas. One such rogue wave caused 
the Ocean Ranger, then the world's largest offshore platform, to capsize in 1982. Rogue 



waves could appear in any place of the world ocean, (see 


Divinsky et aZ.|2004 Earle|1975 


Gutshabash & Lavrenov 


1986 


Irvine & Tilley 


1988 Mori et al. 


2002 


Mallory 


1974) 



They present a threat even to large ships and ocean liners. Until recently rogue waves 
were mostly known to exist only through anecdotal evidence provided by those who had 
encountered them at sea; mariners' testimonies of rogue waves were often treated with 
disbelief despite the damages inflicted by such waves on ships. It all changed after the 
scientific measurement positively confirmed the existence of a large wave at the Draupner 
platform, in the North Sea on January 1, 1995. Since then numerous accounts of rogue 
waves have been reported, mostly in media. It's not just ships and offshore structures 
that need to worry about rogue waves; the U.S. Navy has expressed concern that some 
Coast Guard rescue helicopters lost at sea may have been struck by rogue waves (see 
Freeze|[2006] ). Can rogue wave shoot up so high as to present a danger to helicopters? 



It has been observed that large/rogue waves often appear in sets of more than one, 
the famous sets of threes known as "three sisters" are often quoted. According to some 
accounts three sisters have been observed in Lake Superior. As unbelievable as the stories 
of rogue waves might have been, even less believable are the stories of the phenomenon 
known among Russian sailors as "devyatiy val" (which may be translated as "the ninth 
super- wave" or "the ninth wall of water") and popularized quite unscientifically by Ivan 
Aivazovsky in his famous painting of the same name. According to the folklore, rogue 
waves may appear in groups among which the ninth wave is the strongest and the most 
devastating, hence the name. Ancient Greeks held similar beliefs only it was the tenth 
rather then the ninth wave that was the most powerful. Is the "devyatiy val" just a fiction 
or may there be a grain of truth in it? 

Rogue waves have now been hypothesized as a cause for many unexplained sea losses 
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over the years including the losses over the Bermuda triangle. But can this explanation 
account for the loss of airplanes? Although not very likely, the aforementioned Cost 
Guard report certainly leaves the possibility open. As far-fetched as it might be, it is 
more realistic than the explanations involving "worm-holes", extra-terrestrial life-forms 
or the paranormal. 

There seems to be a fad nowadays to refer to any large oceanic/lake wave as a "rogue 
wave" . Given the abundance of atmospheric, oceanographic and geological phenomena 
on Earth it is very unlikely that all large waves have the same cause(s) and exhibit 
the same behavior. There are currently many models claiming to describe rogue waves, 
at least some must be fairly successful in describing some type(s) of large waves. Since 
large waves considerably differ from each other in shape and behavior it is very unlikely 
that one single model can describe all of them. In this article the term "rogue waves" is 
restricted to the large waves which appear seemingly out of nowhere in relatively calm 
waters, have a relatively short life-span and are localized in a relatively small region of 
space. 

We will show that singular solutions of the Kadomtsev-Petviashvili equation provide 
a surprisingly good description of the shape of at least some large non-rogue waves and 
lead to the prediction of the existence of true rogue waves. We will not address the cause 
of any of these waves. The Kadomtsev-Petviashvili equation is chosen as the simplest 
nonlinear equation of water waves in R 3 many of whose solutions may be computed 
explicitly. One, of course, should understand that the Kadomtsev-Petviashvili equation 
may provide only a limited description of the phenomena. It is quite surprising that the 
Kadomtsev-Petviashvili equation turns out to be so useful as both the singular solutions 
and the large waves exit on the borderline of applicability of the Kadomtsev-Petviashvili 
equation. 



2. Derivation of the Kadomtsev-Petviashvili equation. 

There are several derivations of the Kadomtsev-Petviashvili equation, for the completc- 
ness's sake we shall re-derive it here. Typically the Kadomtsev-Petviashvili equation is 
derived for shallow water waves; however, as the point 4 of the derivation below shows 
the validity of the Kadomtsev-Petviashvili equation may be extended to waves in non- 
shallow stratified water with the vertical component of the motion restricted only to the 
top stratum. To do so consider the motion of fluid, subject to the following assumptions: 

1. The fluid is considered in the region of the three-dimensional space R 3 described by 



where h — constant > 0, and e = constant > is sufficiently small. The set (x',y') G ft, 
z' = —h is the bed of the fluid and the function z' = ehr](x' ,y' ,t') represents the ele- 
vation of the free surface above the reference plane z' = at time t'. We shall refer to 
variables x',y',z',t' as physical coordinates , to ehrj(x / ,y' ,t') as physical elevation and 
to x',y',t',z' = eh-q as physical variables. 

2. The fluid is irrotational, i. e. the velocity vector u'(x' , y' , z', t') satisfies V x u' = for 
all t'. Consequently the fluid has a velocity potential <fr, i.e 



(a/, y') e tt C R 2 , -h< z' < ehn(x' , y' , t'), 



(2.1) 



u' = eV<j>. 



(2.2) 
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3. The fluid is incompressible, i. e. the velocity vector u' and velocity potential <j) satisfy 
' d 2 d> d 2 6 d 2 o 



V • u' = e 



dx' 



dy' 2 dz 



,1 2 



in the region - h < z' < ehrj(x' ,y/,t'). (2.3) 



4. Typically the derivation of the Kadomtsev-Petviashvili equation is based ont he as- 
sumption that the bed is impermeable, i.e the velocity and potential satisfy 



on z' = —h. 



(2.4) 



However, equation (2.4) is also valid if we assume that the flow is horizontal below the 
plane z' — —h, i.e u' ■ n — 0, for z' ^ —h. 

5. On the free surface z' = ehn(x' ,y' ,t') we require that the fluid flows along the free 
surface z' = ehr)(x' ,y' ,t') without ever leaving the surface, i. e. the vertical component 
of the velocity satisfies 

dn d(f> dn d<fi dn 
' dx' dx' 



dz' 



=h 



=h 



dt' 6 dx' dx' ^~ € dy' dy' 



on z'=ehr)(x',y',t'). (2.5) 



6. We assume that the density of the fluid is constant, i. e. 

p = constant. 



(2.6) 



7. The fluid is inviscid, i. e. there is no frictions and the only forces present are gravity 
and pressure. On the free surface z' = ehn(x', y', t') the pressure is given p = p a — 

sV • n = p a — s ( + -^—) where p a is atmospheric pressure, s is surface tension, n 

1 ehn x > x >(l + e 2 h 2 r] 2 ,) l 

is the unit normal to z' = ehrjfx ,y ,t ), — = — ^ and — = 

^1 (l + e 2 h 2 n 2 x ,+e 2 h 2 r 1 2 y ,) 3/2 R * 

ehn v i v >(l + e 2 h 2 rj 2 ,) . , ,.. 

are the reciprocals of the curvature radii. We may assume 



(l + e 2 h 2 V 2 x ,+ e 2 h 2 V 2 y ,) 3/2 

p a = and use assumptions (2.6) to simplify the last two expressions to — - = ehi] x > x ' 
1 

and — — = ehr\yi y i which in turn simplifies the equation for p to 
R2 

p w -seh(ri x > x > + r\ y >y>) at z = ehn. (2.7) 

8. The potential (f> satisfies the Bernoulli's equation e ^7+ e2 ^ |V0| eh(r] x > x > +ri y > y >) + 

d(b 

gz' — C(t'), where C(t') is an arbitrary function. Allowing — to absorb the term 

—C(t) we may rewrite the restriction of the Bernoulli's equation to the free surface 
z' = eht](x', y' , t') as 

S + l 6 ^^ 2 - ^(Vx'x'+Vy'y') +hgr] = on z' = ehr)(x' ,y' ,t'). (2.8) 



9. The fluid motion is essentially one-dimensional with the main motion occurring along 
the x' — axis, the motion and variations along the y' — axis are smaller than the motion 
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and variations along the x' — axis, we symbolically write it as 



d_ 
dy~> 



_d_ 
~dx~' 



(2.9) 



For e sufficiently small the first approximation is obtained by setting e = 0. Equations 
(2.3), (2.4), (2.5), (2.8) become 



V 2 c/> = 
d<fi drj 



-h < z' < r](x',y',t') 
at z' = —h 
at z' = 



90 s . . , 

— + .g?7 - -(r/ x / x / + r^v) = at z = 



(2-10) 
(2.H) 
(2-12) 
(2-13) 



We may look for a solution of system (2.10) - (2.f3) in the form 

V(x',y',t') = e W*'+iv'-»f) } (2.14) 

and 

<j>{x',y',z,t') = f kl (z)e^ kx ' +ly '~^. (2.15) 
Substituting (2.15) into (2.10), we obtain equation — (k 2 + l 2 )fu = whose general 
solution is fki(z) = Cie^ k2+P z + C 2 e~^ k2+l2 z , which upon substitution into (2.15) gives 



us 



<f>(x', y', z, t') = (de^ z + C 2 e-^™ z ) e <(^'+'»'-« *'). 



(2.16) 



Here Ci and C2 are constants of integration which can be determined from boundary 
conditions (2.11) and (2.12) 



Ci = - 



C 2 = - 



-2hVk 2 +V 2 



yjk 2 + l 2 (l-e^vWy Vfc 2 + ^ 2 (l-e- 2 ' l ^+ F ) 

Substituting (2.17) into (2.16) we obtain 

cosh 



(2.17) 



<f>{x',y',z,t') 



(z + h)s/W^7< 



e i(A:x'+Zi/-wt') 



(2.18) 



Vfc 2 + I 2 sinh (Wfc 2 + ? 2 ) 
Boundary condition (2.13) gives us the relationship between w and Vk 2 + I 2 in the form 



to 



2 = Vfc 2 + 1 2 



g + -(k 2 + l 2 ) 
P 



tanh (h^Jk 2 + l 2 ^ (2.19) 



Solving (2.19) for positives and using Taylor approximation along with I <C k which is 
due to (2.9), we obtain 



u= JVk 2 + l 2 



9+-(k 2 + l 2 ) 
P 

\ I 2 h 2 
k+ 2k + Y 



tanh(/i v / fc 2 T] 2 
3s \ 



gph 2 J 



+ lower order terms 



(2.20) 
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As the first approximation to n we may consider 

n = eA^'^'-^+^^^y] *'} + o(fc 3 ) = e **(*'-V9* *') + (jfc), (2.21) 
whose first term satisfies 



a^ 7 



1 gp/i 2 - 3s 

^Vt' + Vx> + 1 7. Vx'x'x' 

gh 6gp 



+ g VVV = (2.22) 



The physical Kadomtsev-Petviashvili equation is obtained by adding to the left hand-side 

dn 2 

a nonlinear term proportional to — — based on some physical argument. The derivation 

ox 

of (2.22) is based on the premise that waves are of the form (2.14)-(2.15) and k, I and 
hy/k 2 + 1 2 are small, i.e. waves (2.14)-(2.15) have long wave length considerably exceeding 
the depth. Due to this the Kadomtsev-Petviashvili equation is often said to describe long 
waves in shallow water. 

Below we discuss a different derivation of the Kadomtsev-Petviashvili equation without 
the requirement that the waves be of form (2.14)-(2.15). Instead we impose conditions 
on the rate of change of the wave profile in time and different spatial directions. We do 
so by introducing dimensionlcss coordinates 

T = <\lr< x " ^r^- Y = %• z = v * - V?' s ~ W S > ,2 - 23 » 

which transform equations (2.3), (2.4), (2.5), (2.8) into 
a 2 $ 9 <9 2 $ <9 2 $ 

6 "dx 2 + W 2 + dz 2 = in region " 1 < z < ^ (2 ' 24) 

^=0 aaZ = -l, (2.25) 

<9<I> dn 9 dn 9 <9<I> dn , 9$ <9n „ 

2 a$ a$ i /a$\ 2 1 /<9$\ 2 i 2 /<9$y 



^-^-5eV,-^+e^ + 2 ^j +2 e^-j +2 ^ ^- j = on Z = en. 

(2.27) 

We look for solutions of system (2.24)-(2.27) in the form 

$ = $o + e$i+e 2 *2 + ---, n = n + en 1 + ... (2.28) 
Substituting (2.28) into (2.24)-(2.25) we obtain 

Qozz + e(*oxx + ^izz) + £ 2 (*im + $oyy + $>2zz) + 0(e 3 ) = 0, (2.29) 
= *iz = $2Z = on Z = -1, (2.30) 

$oz + e(*iz+J7ox) +e 2 ($2Z +mx -Vot -QoxVox') + 0(e 3 ) = 0, on Z = en, (2.31) 
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^$^-$ox+e(j7i-577oxx-$ix + $oT + *oz*iz + ^*ox)+°( e2 ) =0. on Z = en. 

(2.32) 
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Equations (2.28) - (2.30) give us 



$o = A (X,Y,T), $! = --A xx(Z+lf + A 1 (X,Y,T) 

I 1 i 

* 2 = —A 0XXXX (Z + 1) A -- 



(Z + 1) 2 + A 2 (X,Y,T). 



(2.33) 



Af)YY + AixA 

Equating the term of order e° in (2.32) to zero and using <j> z = A z = we obtain 

Vo = $ox = A ax . (2.34) 

Equating the coefficient of e in (2.32) to zero and using $ z = A^z = 0, Z = 1 + en as 
well as (2.34) we obtain 

Vi-Svoxx-^ix+^ot + ^Ix = m-Smxx-A-Lx + ^mxx+AoT+^Vo = °- ( 2 - 35 ) 

Substituting $ z = A oz = 0, Z = 1 + en, $ lz = -A oxx , $ 2 z = \.A 0XXXX - A 0Y y - 

6 

A\xx into equation (2.32) we obtain that the coefficient of e° is trivially zero while 
equating the coefficient of e to zero and using (2.34) we obtain 

-AqxxVo + $2Z+Vix - Vot - ^oxVox = 

1 (2.36) 

- ZiloxVo + j,mxxx - Aqyy - Aixx + Vix - Vot = 0. 
6 

The differentiation of (2.35) with the help of 770 = A ox further gives us 



Vix - Sr] 0XXX - Ai XX + -Voxxx + Vqt + VoVox = 0. 



Subtracting (2.36) from (2.37) we obtain 

1 



2r?oT + 



S Voxxx + 3r] mx + A aY Y = 0. 



Differentiating (2.38) and using (2.34) we obtain 



d_ 
dX 



2rj, 



0T 



g - S ) Voxxx + 3voVox 



V0YY = 0, 



or in dimensional coordinates 
d 



dx' 



1 pgh 2 — 3s 3e 

Vot< + Vox' H -p, Vox'x'x' + -prVaVox 



1 



:V0y'y> 



0. 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



fgh "" 6pg 

which we shall call physical Kadomtsev-Petviashvili equation. Notice that even if we 
dropped n Y Y in (2-27), or, equivalently, assumed that R 2 = +00 in the derivation of 
(2.7), r/oYY would still appear in (2.39). 
The change of variables 

4/ 



V 



o 2pg . 3a 2 t . 2 , 



or equivalently, 



a 2 = 



2pg 



t = 



pgh 2 -3s' 3a 2 

turns equation (2.40) into 

d_ 
dx 



gh t' . i—r . , r- . 3ea 2 

x = x - \Jgh t, y = y V 2 , / = — — 770, 



ft + fxxx + 6.ffx + 3a 2 f yy = 0, a 2 =±l, 



(2.41a) 



(2.41b) 



(2.42) 
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which is known as the Kadomtsev-Petviashvili equation or simply KP. Notice that the 
physical Kadomtsev-Petviashvili equation given by (2.40) differs from the Kadomtsev- 
Petviashvili equation given by (2.42). It is the physical Kadomtsev-Petviashvili equation 
(2.40) that describes the physical elevation ehn(x' , y' , t') in terms of physical coordinates 
x',y',t'; equation (2.42) describes a nonphysical quantity f(x 7 y,t) in terms of variables 
x, y, t in the frame of reference moving with velocity \/g~h in the physical frame of reference 
x',y',t'. 

Notice that formulas (2.2) and (2.34) imply that for z' > — h 

u' x , = e— = ey/gJ <f> x + o(e) = e^/g~h ri + o(e), (2.43) 

where, according to (2.2), u' x , is just the x' — component of the velocity of the flow. Thus 
the physical Kadomtsev-Petviashvili equation may be viewed as an equation not only for 
elevation 770 but also for the x'— component of the velocity of the flow u' x , . Moreover, the 
first two formulas in (2.33) and formula (2.34) imply that for z' > —h 

u' z , = eg = Jgh? $ lz +o ( e 3 / 2 ) = mx +{e^) = -e{h+z')^h ^+0 (e 3 / 2 ) , 

(2.44) 

d 

keeping in mind that — - = 0(y/e) due to (2.23). 

Formulas (2.43), (2.44) give us the components of the velocity of motion of the fluid 
not of the profile ehrj(x / ,y' ,£'). Whereas the profile ehn may move in one direction the 
fluid on the surface z' = ehn(x' , y' , t') may move in another direction. 



3. Special solutions of KP. 

The simplest solution of KP is a soliton wall solution 

»2 / Jp+q)x + (q 2 -p 2 )y-(p+q)[(p+q) 2 +3a 2 (p-q) 2 ]t\ 

fix, y,t) = 2- 2 in (^1 + — j , (3.1) 

where c, p, q are real constants. Nonlinear superposition of TV soliton walls is given by 

d 2 

f(x,y,t) = 2— IndetA (3.2a) 

where 

c e (Pn+q„)x+(ql-p 2 n )y-(p„+q n ) [(p n +q n ) 2 +3a 2 (p n -q n ) 2 ] t 

Ann = 5 mn + — , m,n=l...N, 

Pn + Qrn 

(3.2b) 

S mn are the regular Kronecker symbols, c„, p n , q n p are real constants. 

We may determine the velocity v = (v x ,v y ) of the motion of a single soliton wall by 

' ; [ (p + q)x + {q 2 -p 2 )y-(p + q)[(p + q) 2 + 3a 2 (p-q) 2 ]t = 
(p + q)v x + (q 2 - p 2 )v y - 4(p 3 + q 3 ) = 0, which gives us 

> V * + 1 , , o 9/ ^ V V = L ( 3 ' 3 ) 



setting df 



(p + q) 2 +3a 2 (p-q) 2 x (p + q) 2 + 3a 2 (p - q) 2 v 

So velocity is not defined uniquely but up to an additive term proportional to the unit 
p-q 1 



vector 



^l + ip-q) 2 ' ^l + (p-q) 2 
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The components v — (v x ,v y ) of the velocity of motion of the superposition of two 
soliton walls satisfy the system of equations 



v x + (qi - PijVy = (pi + qi) 2 + 3a 2 (pi - qi) 2 

V x + (<?2 ~ P2)Vy = (P2 + <?2) 2 + 3a 2 (p 2 - <?2) 5 



(3.4) 



and are uniquely defined provided the determinant of the system is nonzero, i.e. p± — qi =/= 
P2 — qi ■ We cannot attach a meaningful definition of velocity to the motion of nonlinear 
superposition of three or more soliton walls. The velocity of the soliton wall v = (v x , v y ) 
is not to be confused with the velocity of the flow given by the vector v! , the former is 
the velocity of motion of the profile while the latter is the velocity of motion of the fluid. 

One may notice that while the derivation (2.14)-(2.22) of the Kadomtsev-Petviashvili 
equation starts with the assumption that solutions are of the form (2.14)-(2.15), functions 
(2.14)-(2.15) themselves are not solutions of KP (see Kovalyov & Bica 2005). In search 
for solutions of KP of the form close to (2.14)-(2.15) they arrived at 



d 2 

f{x lV ,t) = 2—\n 



T = p + x cos(ax) + 2 



2AT - sin 2r 
A sin(ax) 



p cos(ax) 



y+ 



12 



A 2 cos(ax) — o?\j? cos(ax) + 2aA^ sin(ax) 



T = 7 + A.t - 2A/iy + 4A(A 2 - 3a 2 p 2 )t, 



(3.5a) 

(3.5b) 
(3.5c) 



where A, /i, 7, p are some constants. If the constants A, /1, x, 7, p are real so is the solu- 
tion (3.5). Nonlinear superposition of such solutions is of the form 



u(x,y,t) 

where K is an N x N matrix with the entries 

/ K u K 12 

K = 



d 2 

2— lndetK 
ox z 



K 



21 



K- 



22 



\ K N i K N2 



K lN \ 
Knn J 



K — T — 

J1 ™ — J- n 



sin 2r ? , 
2A n 



(3.6a) 



(3.6b) 



(3.6c) 



(An ~ Afc) sin(r„ - T fc ) 
a 2 (p n - Hk) 2 + (A„ - Afc) 2 

(Pn - Hk)cOs(T n + Tfc) 

a 2 (/i„ - p k ) 2 + (A„ + A fc ) 2 
A„ sin(axn) 



Tf„ = p n +x cos(ax„) + 2 



a 



(A n + Afc)sin(r ra + r fc ) 
a 2 (/z„ - p k ) 2 + (A„ + Afc) 2 

(M« - Pk)cos(T n - T k ) 
2 (^„-Mfc) 2 + (A„-A fe ) 2 



p n cos(axn) 



n ^ fc(3.6d) 



12 A 2 cos(ax„) - a 2 pl cos(ax„) + 2aX n p n sin(ax„) t, (3.6e) 

T„ = 7„ + X n x - 2X n p n y + 4A„(A 2 - "ia 2 p 2 n )t, (3.6f) 
where A„'s,/i„'s,Xn's, 7n's, p„'s, are constants. Change of constants p n — > — p n , \n 
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Figure 1. Graph of regularization (3.7) of (3.5) with a=l, A = 0.5, fi=—0.1, 7 = 0, p=0, 

x=o, t=o. 



h 7r allows us to change the sign in front of T„ while the change 7„ — ¥ 7„ ± — for all 

a 2 



(A n + A fe ) sin(r„ + r fc ) 



(fin - fik) cos(r„ + r fc ) 



(fJm - Vk) 2 + (A„ + \ k ) 2 ' 



n = 1, 2, 3 . . . , N allows us to change the signs in front of all , 

or 

a 2 (ii n - Hk) 2 + (A„ + A fc ) 2 ' 

Solution (3.5) is singular and approaches — oo as we move towards the singular curve 
2AT — sin 2r = 0. To be physically plausible it needs to be regularized to remove infinities; 
since there is no clear way to determine what kind of regularization is most suitable we 
choose an ad-hoc regularization 



f -> F 



J_ 
l/l 



In 



In 



due to its properties 

if0<|/|<l; 



( e f - 1 + l) + 



0, if / = 0; 



F4 - ln(ln 2 + 1), as / -> -oo; F In / as / 



(3.7a) 

(3.7b) 
(3.7c) 



Conditions (3.7b) assure that the regularization is close to the original function in the 
domain of validity of KP , whereas conditions (3.7c) remove negative infinity and reduce 
large amplitudes to make them more physically plausible. 

Regularization (3.7) of (3.5) with \— is shown in Figure [l] while regularization (3.7) 
of (3.5) with x + 1 is shown in Figure [5] 

We may determine the velocity v = (v x ,v y ) of motion of the profile given by (3.5) by 

(1 f Aa; - 2A ) uy + 4A(A 2 - 3a 2 ^ 2 )tj = \v x ~ 2Xfxv y + 4A(A 2 - 3a 2 (j?) = and 



setting 
d " 



df 



dt 

cos(ax)x+2 



cos(a\)v x +2 



Asin(ax) 
a 

A sin(ax) 



-/j cos(a%) 
-[i cos(ax) 



y +12 
;„ +12 



0, which give us a system of equations 



A cos(ax)— ol jj, cos(ax)+2aXfi sin(ax) 



A cos(ax) - a f 1 cos(ax)+2aA/xsin(ax) = 



A tan(ax) 



/' 



12 



A — a /i + 2aA/i tan(ax) 



(3.8a) 



v x - 2[iv y + 4(A 2 - 3a 2 (i 2 ) = 0. 





If A tan a% 7^ the determinant of system (3.8a) is nonzero and the system has a unique 
solution 



Ua- = —12a fi — 4A , w.„ = —12a /1 . (3.8b) 

tan ax tan ax 

Far away from the singular curve the oscillatory portion of the waves moves in the direc- 

.. ( 1 -2/i \ , ( v x ~2fiv y 12a 2 n 2 -4\ 2 
tion — i — with the speed of — = ^ — which is 



vV /TTv' yiTv/ ' ' y/1 + V V 1 + V 

independent of x- The singular curve itself stays within distance — , -p from the 

2|Acos(ax)| 

straight line P T — r +x + 2 [ - tan ^ ax ^ - J y+ 12 ["a 2 -a 2 u 2 + 2aAu tan(ax) 

cos(ax) 

moves in the direction of 



1 2 



12 

speed of 



A 2 - a 2 /^ 2 + 2aA/^tan(ax) 




1 + 4( M - 



A tan ax ^ 



bounded for all values of x- For x = profile 



(3.5) moves in the direction of ( — , — = | with the speed of — (^_i^__J. 

\yi + V v / i+¥/ \/i + V 

while at the same time oscillating along the same direction. The time evolution of such 
a wave is shown in Figure [3j it is similar for a = 1 and a = i. 

We cannot attach a meaningful definition of velocity to the motion of nonlinear super- 
position of two or more solutions (3.5). The velocity v = (v x ,v y ) is not to be confused 
with the velocity of the flow given by the vector u', the former is the velocity of motion 
of the profile while the latter is the velocity of motion of the fluid. 

Each such solution consists of two components separated by a singular line. We may 
retain the component on one side of the singular line and set the function to zero on the 
other side 

f(x, y, t) given by (3.5a), if 2AT - sin2r > 0, 
f(x,y,t) = { (3.9) 
0, if 2AT-sin2r < 0. 
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Figure 3. Graph of regularization (3.7) of (3.5) with a = l, A = 0.5, n=— 0.1, 7 = 0, p = 0, 
X = 0, —30 < x < 15, —30 < y < 30, for the values of t shown. 



The obtained function is singular along the line 2AT — sin 2r = and thus could not 



be expected to model any physical wave in a neighborhood 



2AT - sin2r < 



< e, e is 



some constant, of the singular line. But away from the singular line in the domain 



2AT- 



sin2r < > e function f(x,y,t) provides a quite good description of undular bores or 
waves with the advancing strongly-pronounced front followed by a train of well-defined 




Figure 4. Graph of regularization (3.7) of a simply connected component of (3.5) with 
a = l, A = 0.5, /x = -0.1, 7 = 0, p = Q, x = 0, t = 0. 




Figure 5. Garonne river tidal bor e at Arcins on July 6, 2008 around 07:10 (views from right 
bank), described at ( Chanson]|2011 1. Courtesy of Hubert Chanson. 



free-surface undulations. Near the singular line f(x,y,t) has a large negative trough 
followed by a large positive crest as shown in Figure [4] Formula (2.43) for the physical 
velocity u' x , — t\fgh t/q + o(e) shows that while water at the crest of the wave continues to 
move forward, water in the trough moves backwards leading to wave overturning which 
cannot be described within the framework of KP. Numerous discussions and pictures 



of undular bores may be found on the Internet, e. g. (Chanson 2010), one of which is 
reproduced in Figure [5] 

An undular bore-like structure was also observed in December 26, 2004 tsunami. As 
the tsunami passed around the Phi Phi Islands and moved towards Koh Jum island off 
the cost of Thailand it took undular borc-likc form and was photographed by Anders 
Grawin. The pictures are reproduced in ( Constantin||2008| and (Anders 2004) with one 
of them shown in Figure [6j 

Several undular bore- like waves are shown in Figure [7] in the image captured by ESA. 

The very existence of the singular line indicates that most of the power of an undular 
bore is concentrated in the trough at the leading edge rather than at the bore's crests. 



The presence of large forces at the leading edge of undular bores is described in ( Simpson 
et q/.|2004 Mouaze et a/.|2010 Wolanski et aL|2004| . The account of an undular bore in 
^ Wolanski et aZ.|2004 ) as " unsteady motion ... sufficiently energetic to topple moorings 
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Figure 6. The tsunami of December 26, 2004 approaching Koh Jum island, off the coast of 
Thailand, after it passed around the Phi Phi islands. Copyright Anders Grawin, 2006. Repro- 
duced from http://www.kohjumonline.com/ anders.htmll with permission. 



mm 
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Figure 7. Image captured by ESA, courtesy of ESA and Deutsches Zentrum fur Luft- und 
Raumfahrt (DLR). On the left is the original image, on the right is the same image with some 
waves of the type shown in Figures [I] and [4] enclosed in black rectangles. The middle of the 
horizontal white line marks the wave with the highest crest, it is the leading crest of a wave of 
the type shown in Figure [4] 



that had survived much higher, quasi-steady currents of 1.8 m/s " also supports the claim 
that the real power of an undular bore is in its leading trough rather than relatively 
modest crests even though the leading trough may not be well-pronounced. As follows 
from (3.5) the leading crest may change its shape. 

Often the undular bore-like waves are misinterpreted as trains of solitons with dimin- 
ishing amplitudes. If that were true, the solitons would be dispersing as the smaller ones 
travel slower than the taller ones. Yet the undular bore-like structures do not disperse 
and preserve their shapes for quite a while behaving more like solutions (3.5) than trains 
of solitons. We refer to solutions (3.5) as "harmonic breathers". 

Although most undular bores-like waves have the shape shown in Figure [4] with the 
strongest crest-trough pair leading the wave, there might have been observations of waves 
with the shape resembling Figure [I] with the strongest crest-trough pair being in the 
middle rather than in the lead of the wave. According to (Dudley 1998), the Aleutian 



tsunami of 1946 produced multi-crest waves with the third or fourth crests being the 
highest and most violent, and a bore at the Waimea Riverr on Kaua'i is reported to have 
the sixth crest as the highest and most destructive. 

Function (3.5) is singular along the curve 2AT — sin2r = and thus could not be 



2AT - sin 2r < 



< e, e is 



expected to model any physical wave in a neighborhood 
some constant, of the singular curve. However, away from the singular curve the function 
describes the waves observed in the river Seven Ghosts in Sumatra, Indonesia and shown 
in Figures [8j [9] taken from the video ( |Rip Curl fmT^ 

Formula (2.43) for the physical velocity u' x , = e^/gh no + o(e) shows that while water 
in the trough moves backwards, at the points behind the singular curve where the water 
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Figure 8. Time frame t = 1 : 21 from the video at ( |Rip Curl~||2011[ ). 





Figure 9. Portion of time frame t = 2 : 11 from the video at ( Rip Curl |2011 1 



level is relatively high water continues to move forward, leading to wave overturning but 
only at some points along the singular curve in agreement with the video taken at the 
river Seven Ghosts. 

Another class of explicit solutions of KP may be obtained by simple formal substitu- 
tions 

A — > iX, x ~~ * iXi 7 *7 m to (3.5) 



A„ -> iA„, Xn -> iXn, In -> «7« into (3.6) 
Substitutions (3.10) give us 



f(x,y,t) = 2-jLhi 



T = p + a; cosh(ax) — 2 



2AT - smb. 2r 
A sinh(ax) 



+ iicosh(ax) 



(3.10) 



(3.11a) 



12 



A cosh(ax) + a fi cosh(ax) + 2aA/isinh(ax) 



t, 



(3.11b) 



r = 7 + Ax - 2X)iy - 4A(A 2 + 3a 2 /j, 2 )t, 



(3.11c) 



where A, /i, x, 7, p are some constants. If the constants A, /x, %, 7, p are real so is the solu- 
tion (3.5). Nonlinear superposition of such solutions is of the form 



d 2 

u(x,y,t) = 2—— lndet K 
ox z 



where K is an TV x N matrix with the entries 



K = 



( K tl K 12 
K21 K 2 2 



K 1N \ 

K 2 n 



(3.12a) 



(3.12b) 



\ K N i K N2 ■ . ■ K NN J 
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Figure 10. Graph of regularization (3.7) of (3.11) with a = l, A = l, /i = 0.05, 7 = 0, p = 0, 
X = 0.6, t = Q. 



T 



K 



nk 



sinh 2r„ 
2A„ 

(A„ - A fe )sinh(r„ - Y k ) (A„ + A fc ) sinh(r„ + Y k ) 



(3.12c) 



a 2 {p n - /i fe ) 2 - (A n - Afc) 2 



2 (Ai«-<Ufc) 2 -(A„ + A fc ) 2 



a 



(Mn - Mfe) cosh(r„ + r fe ) (fi n - /i fe ) cosh(r r , 



r,) 



a 2 (p n - p k ) 2 - (A„ + A fe ) 2 a 2 (p n - p k ) 2 - (A„ - Afe) 2 _ 
A n sinh(axn) 



p n + xcosh(axn) 



p n cosh (ax,, 



12 



A 2 cosh(axn) + a 2 fi 2 cosh(axn) + 2aA„^„ sinh(axn) 



t. 



n^k (3.12d) 



(3.12e) 



r„ - 7« + Kx - 2X n p n y - 4A„(A 2 + 3a 2 ^ 2 )i, (3.12f) 
where X n 's,fjL n 's,Xn's, 7n's, p n 's, are constants. Change ol constants p n — > — p n , \n — > 

— +in allows us to change the sign in front of T„ while the change 7„ — > 7„ ±i— for all 
a 2 

i o q M I, . -i . , . . , , t (An + A fc ) sinh(r„ + r fc ) 
n = 1, 2, 3 . . . , Jy allows us to change the signs in front ol all — — — -z , 

er(/U„ - Hk) - (K + Afe) 2 

{pn - Mfc) cosh(r„ + T k ) 



HPn - Mfe) 2 - (K + -^fc) 2 ' 



Due to (2.43) the fluid moves to the right whenever / > 0, such regions are shown in 
Figures [lOj [IT] in light shades; the fluid moves to the left whenever / < 0, such regions 



are shown in Figures 10 11 in dark shades. When for x 7^ the fluid moving to the 
right collides with the fluid moving to the left, we obtain a "crossing" of light and dark 
lines. The velocity v — (v x , v y ) of the motion of the "crossing" is the same as the velocity 
of the wave profile and its components may be uniquely determined from the system of 
equations 



v x - 2 



A tanh a\ 



Vy = 12 



A 2 + a 2 p 2 + 2Xpa tanh a\ 



(3.13a) 



2pv v = 4(A 2 + 3a> 2 ) 



obtained from (3.8a) by means of substitutions (3.7). If the determinant of system (3.8a) 
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Figure 11. Graph of regularization (3.7) of (3.11) with a = i, A = 0.6, /i = 0.05, 7 = 0, p = 0, 
X=0.6,t = 0. 



is nonzero, i. e. Atanhax ^ 0, the system has a unique solution 

v x = 4A — 12a fj, , v v = 12a /i. (3.13b) 

tanh a\ tanh ax 

m , , v x - 2/j,v v 12a 2 /i 2 + 4A 2 . . , , . 

The velocity component — = ^ — is independent of x while the 

component — x v becomes infinite if x — 0. F° r X 7^ t ne motion of profile 

(3.11) may be viewed as superposition of two motions: one is the transverse motion 

■ , ■ r / 1 2/1 \ . , , , r 12c*V +4A 2 , , 
m the direction of — - , with the speed oi , — and the 

other one is the lateral motion of the "crossing" along the wave in the direction of 

— = , — -j^=^= | with velocity ^ Vx = . For x = profile (3.11) moves in 

v / i + V \A+ V/ V 1 + V 

1 2u \ . , , , r 12a 2 Ai 2 +4A 2 
the direction oi — . , . with the speed of . — while at the 

same time the crest grows, jumps over the singular line and then disappears, in this case 
the velocity of the lateral component of the motion becomes infinite. The time evolution 
of such waves for a = 1 and a = i is shown in Figures [12] and |13| 

We cannot attach a meaningful definition of velocity to the motion of nonlinear super- 
position of two or more solutions (3.11). The velocity v = (v x ,v y ) is not to be confused 
with the velocity of the flow given by the vector u', the former is the velocity of motion 
of the profile while the latter is the velocity of motion of the fluid. Each of the Figures 
[TUJ [TT] consists of two long crests preceded or followed by long troughs, each crest-trough 
pair may be viewed as a separate wave looking somewhat like what is shown in Figures 

Due to formula (2.7) the pressure at the point in the "crossing" given by p+x cosh(ax)— 
" A sinh(ax) 



/tcosh(ax) y — 12 A cosh(ax) + a M cosh(ax) + 2aA/i sinh(ax) 



t 



a 

7 + \x — 2\py — 4A(A 2 + 3a 2 /i 2 )i = is +oo, thus the pressure in the physical wave 
modeled by (3.11) should be much larger at the "crossing" than anywhere else. 
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Figure 12. Graph of rcgularization (3.7) of (3.11) with a = l, 
X = 0, —20 < x < 20, —35 < y < 35, for the values of t shown. 



A = l, 



-0.05, 7 = 0,p = 0, 



The physical phenomena described by such functions are waves that stretch along an 
almost straight half-line, with amplitude sharply dropping at the end point and slowly 
decaying to zero in the other direction, one such wave shown in Figure |16| Much better 



illustrations may be found at numerous web sites, e.g. (see Videos||2010[ |Riptide 2010). 



An image obtained by ESA from a satellite and reproduced in Figure 17 exhibits several 
such waves, each is shown by a pair of white and black strips next to each other. One 
such wave is described in ( Grahairi||2000 1 as " Only the very long swell, of about 15 feet 
high and probably 600 to 1000 feet long. . . . We were on the wing of the bridge, with a 
height of eye of 56 feet, and this wave broke over our heads. ... we were diving down off 
the face of the second of a set of three waves, so the ship just kept falling into the trough, 
which just kept opening up under us." 

Although some of the literature refer to such waves as "rogue waves", we refrain from 
using the term for this type of waves as typically such waves appear in "predictable" 
situations unlike real rogue waves that appear seemingly out of nowhere. While the 
sharp crest of the wave is clearly visible the trough may not be as apparent. 



18 



M. Kovalyov 



t=-1.8 


t=-0.04 

/ 


f 


I 




I 




t=0.04 







-0.5 0.0 1.0 2.0 

Figure 13. Graph of regularization (3.7) of (3.11) with a = i, A=l, fJ.= — 0.05, 7 = 0, p=0, 
X = 0, —20 < x < 20, —35 < y < 35, for the values of t shown. 




-20 



Figure 14. Graph of regularization (3.7) of the upper half of (3.11) with a = l, X=l, fi — 0.05, 
7 = 0, p=0, X = 0.6, t=Q. 



Another type of a large wave may be obtained by substituting 



7T 7T 

a = l, p^f-ip, x^x + i-, 7^7 + *^; 



(3.14) 
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Figure 15. Graph of rcgularization (3.7) of the upper half of (3.11) with a = i, A = 0.6, (i = 0.05, 
7 = 0, p=0, X = 0.6, t=Q. 




Figure 16. Physical wave with the shape predicted by Figures [T4]-|15| 
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Figure 17. Image captured by ESA, courtesy of ESA and Deutsches Zentrum fur Luft- und 
Raumfahrt (DLR). On the left is the original image, on the right is the same image with some 
waves of the type shown in Figure [14] enclosed in black rectangles and two waves of the type 
shown in Figure [10| enclosed in white rectangles. In the latter the "crossings" can be clearly 



into (3.11), which gives us 



2AT + cosh2r 



T=p+x sinh \— 2 A cosh x+M sinh \ V — 12 A 2 sinli x+M 2 sinh x+2A/i cosh \ t, 



(3.15) 



r = 7 + Xx - 2Xny - 4A(A 2 + 3fi 2 )t. 
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Figure 18. Graph of regularization (3.7) of (3.15) with a = l, A = 0.4, ^ = 0.05, 7 = 0, p = 0, 
X = 0.6, t=Q. 



71 . 



where x, A,/i,/9 are rea l constants, 7 is either a real constant or a real constant i^*- 



Regularization (3.7) of (3.15) is shown in Figure 18 



Due to (2.43) the fluid moves to the right whenever / > 0, such regions are shown 
in Fig. 2 in light shades; the fluid moves to the left whenever / < 0, such regions are 
shown in Figures [TO] - [IT] in dark shades. When the fluid moving to the right collides with 
the fluid moving to the left, we obtain a "collision" of light and dark lines. The velocity 
v = (v x , v y ) of motion of the "collision" is the same as the velocity of the motion of the 
wave profile, its components may be uniquely determined from the system of equations 



v, - 2 



A coth x + A 4 



Vy = 12 



A 2 + n A + 2A^tcothx 



(3.16a) 



v x - 2fiv y = 4(A 2 + 3^ 2 ) 



obtained from (3.13a) by means of substitutions (3.14). Since the determinant of system 
(3.16a) is always nonzero, the system has a unique solution 



4A 2 - 12^t 2 - 8A^tanhx, 



-4Atanhx - 12/i. 



(3.16b) 



The velocity is finite for all values of x- We cannot attach a meaningful definition of 
velocity to the motion of nonlinear superposition of two or more solutions (3.15). The 
velocity v = (v x , v y ) is not to be confused with the velocity of the flow given by the vector 
u 1 ', the former is the velocity of motion of the profile while the latter is the velocity of 
motion of the fluid. 

Waves of type (3.15) would be extremely unstable, so it would be very unlikely to see 
one as shown in Figure [18} However the portions shown in Figure [19] should be observable 
at least for short periods. Indeed, Figures [20[[22] show such waves. 

Formula (2.43) implies that the higher up are the points of the fluid the faster they 
move in the x direction ultimately leading to wave overturning which destroys the wave. 
However, neither solutions (3.11) nor solutions (3.15), (3.18) exhibit wave overturning 
due to the limitations of KP equations. This limits the applicability and usefulness of 
solutions (3.11)-(3.15) yet they do describe certain aspects of the corresponding physical 
waves such as sharp crests, the presence of troughs, the motion as a whole, etc. The 
life-span of the corresponding physical waves might be relatively long or fairly short 
depending on how long conditions (2.1)-(2.9) are sustained. 
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Figure 19. Central portion of the wave shown in Figure [18] with terminology used in the text. 




Figure 20. Physical wave with the shape predicted by Figure [19] U-shaped crest, depression 
inside the crest and U-shaped trough are well-pronounced; bulging inside the trough is seen but 
not well-pronounced. The picture is a frame from a video shot near Kiama on the NSW South 
Coast. Appeared in ( Ackerman|2011 1. 



The scaling transformation 



x^xS, y^yS 2 , t^td 3 , f -> (3.17b) 

o 6 

A->^, (3.17a) 

transforms each of the solutions of KP considered above into a solution of KP of the 
same type, the velocity undergoes transformation 

V X Vy 

v x -> , v y ^ (3.17c) 
Transformation (3.17) allows us to re-scale the graph of the corresponding solution. 



equivalent to 
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Figure 21. Physical wave with the shape predicted by Figure [19] U-shaped trough and bulging 
inside the trough are seen; U-shaped crest and depression inside the crest are barely seen. The 
original photo was taken by Seamus Makim near El Fronton, Canary Islands in 2010. 




Figure 22. Physical wave with the shape predicted by Figure [19] U-shaped trough and bulging 
inside the trough are well-pronounced; U-shaped crest and depression inside the crest are not seen 
and may nor may not be present . The original photo was taken by Mike Maxted in south-west 
Western Australia. Appeared in ([Maxted 20101. 
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4. The true rogue waves. 

In the previous section we discussed the use of some singular solutions of KP to model 
large waves. They may also be used to provide a possible explanation for true rogue waves. 
Indeed, consider superposition of two waves of type (3.11) given by formulas (3.12). Each 
wave moves in a corresponding direction with its " crossing" traveling laterally along the 
wave as discussed after formula (3.13b). When the "crossings" of the two waves come 
close to each other they produce bi-lateral compression which literally squeezes water to 
the top by the force exhorted by the "crossings"; a very similar effect is obtained if one 
claps hands under water just below the surface. This results in a drastic short-term jump 
in amplitude localized in a small region much the same way true rogue waves do. We refer 
to such a jump as OTIN which is an acronym for One Time INtense big wave/splash/wall 
of water/etc. Bi-lateral compression may create a very large OTIN even if the two waves 
themselves have rather modest amplitudes, making it appear as if the OTIN rises out of 
nowhere. 

To remove from plots physically unattainable large positive and negative values from 
f(x,y,t) given by (3.12) let us use ad-hoc re-normalization 



r / 



in 



In 



f-1 



f M (x,y,t) = < 



l/l 

f, if < / < M, 
M, if / > M, 



1+1 



if / < 0, 



(4.1) 



where M is a sufficiently large positive number. Re-normalization (4.1) completely cuts off 
the values of f(x, y, t) above M and provides a better resolution in the range < / < M. 
The previously-defined re-normalization (3.7) does not cut off any values and allows us to 
look at all of f(x, y, t) > by using a logarithmic scale yet the very use of the logarithmic 
scale reduces the resolution. 



Figure 23 shows re-normalization (4.1) with M = 10 of the evolution of of wave (3.12) 
with N = 2, a=l, xi =0.6, A x =0.5, ^=0.2, 7l =0, Pl =0, xa=-0.7, A 2 =l, ^2=0.5, l2 = 
0, /92 = in the region —15 < x, y < 15 . The height of the surface stays relatively low 
for t < -0.23 and t > 0.23 while at the time interval -0.23 < t < 0.23 a powerful OTIN 
rises from the water seemingly out of nowhere only to disappears seemingly into oblivion. 

Under certain conditions an OTIN may split into several parts. To demonstrate it 
consider wave (3.12) with N — 2, a = l, xi = 0.6, Ai=0.5, /ii=0.01, 7i=0, pi=0, \2 = 
—0.7, A2 = l, /it2=0.5, 72=0, P2=0, which diffcrers from the previous example solely by 
the value of /zi. Time evolution of re-normalization (4.1) of this function is shown in 
Figure [24] in a moving coordinate frame. Notice that the OTIN appears to be the largest 
and the most powerful in the beginning before it breaks up and in the end after all the 
parts are recombined together. The event may appear to the crew of a vessel carried by 
the current as a very powerful wave followed by a few less powerful ones with another 
powerful wave int he end. Whereas the first one hits the vessel which is still intact and 
with the relatively fresh crew it may not appear that strong, the last one hits the vessel 
already heavily pounded with an exhausted crew. It is the last one that would appear to 
be the most devastating to the crew thus leading to the legends of the " devyatiy val" . 



Figure 24 also shows that the small white spots indicating high elevation cluster in 
groups of threes, e. g. the lower three white spots in the frame t = —4 , the upper three 
white spots in the frame t = 4, the upper three white spots the lower three white spots 
in the frame t — 0. Some rogues waves have been reported to appear in groups of threes, 
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Figure 23. Regularization (4.1) with M = l of time evolution of (3.12) with N = 2, a=l, 
Xi =0.6, Ai =0.5, /ii =0.2, 71 =0, pi =0, X2 = —0.7, A2 = 1, [1.2 = 0.5, 72 = 0, P2 = in the region 
-15 < x,y < 15. 



e. g. on April 16, 2005 a cruise liner Norwegian Dawn encountered a series of three 21.34 
m rogue waves off the coast of Georgia; in March 2012 a cruise liner Louis Majesty 
encountered a series of three 7.9 m rogue waves between Cartagena and Marseille; there 
are reports of rogues waves appearing in groups of three on Lake Superior. 



OTINs rise in the regions shown in Figures 23 - 24 in black or close to black shades, 



these are the regions whose surface is considerably below the sea level and where a lot 
of underwater traffic takes place, such regions are just giant troughs in the ocean. Since 
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Figure 24. Regularization (4.1) with M = 1 of time evolution of (3.12) with N 
Xi=0.6, Ai=0.5, mi=0.01, 71=0, pi=0, x 2 =-0.7, A 2 =l, M2=0.5, 72 =0, p 2 =0 



2, a=l. 



the surface of such a trough is below the sea level, it makes even a modestly tall OTIN 
to appear considerably taller, e.g. an 4 meter OTIN rising from a trough 3 meters deep 
appears to the crew of a ship that happens to be in the same trough to be 7 meters tall 
thus contributing to more drastic perception of the OTIN. 

Both graphs have portions somewhat resembling linear waves. Figure [25] provides con- 
tour maps of the portions < / < 0.6 of two frames of Figure [23] (in the first row) and 
two frames of Figure 24 (in the second row). These waves appear before, after and around 
the OTINs contributing to the appearance of calm waters and exacerbate the impression 
of OTINs appearing seemingly out of nowhere. 



5. Conclusion and final remarks. 

Large ocean waves have always been the subject of human fascination and fear. Yet 
the shear nature of such waves makes it extremely difficult to created a good mathe- 
matical model. In §3 singular solutions of KP are used to describe some aspects of some 
large waves; the description is far from being complete yet to the best of the author's 
knowledge there is no better one. Such limited models are often referred to as "mathe- 
matical caricatures" due to the fact that just like regular caricatures they provide rather 
sketchy description of the phenomena with some features exaggerated and some features 
suppressed. It is rather remarkable that the singular solutions of KP actually do describe 




some properties of large waves as both the singular solutions and large waves are near 
the borderline of applicability of KP. In §3 singular solutions and their superpositions 
are given by formulas (3.1), (3.2), (3.5), (3.6), (3.11), (3.12), (3.15), yet so far the author 
has not succeeded in combining these formulas into a single one. 

In §4 we discussed the phenomenon of OTINs that seems to provide an explanation 
of rogue waves that suddenly appear seemingly out of nowhere and after a short while 
disappear without a trace. OTINs also seem to explain why rogue waves are often encoun- 
tered in clusters of three. Formulas and illustrations of §4 suggest that OTINs appear due 
to the bi-lateral compression of fluid under the ocean surface caused by two converging 
"crossings". It also provides a method of predicting where OTINs may rise. To determine 
the danger zones one needs to look at the regions of intersections of giant troughs travel- 
ing across oceans, such regions should be detectable from satellites. The description and 
explanation of rogue waves provided in §4 considerably differ from other existing models. 

Although the mathematical formulas discussed here provide rather limited description 
of the corresponding physical waves, they do so in a coherent manner providing a unified 
albeit very primitive theory of large waves predicting the existence of "crossings" for the 
waves shown in Figure |10| the existence of bulge inside the trough for the waves shown 
in Figure 19 the existence of OTINs, the clustering of rogues waves in groups of three 
as shown in Figure [24} etc. In other words, they show that there "a method in the mess" 
and there are underlying principles responsible for all these phenomena. 
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